Concerning the quark condensate 
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A continuum expression for the trace of the massive dressed-quark propagator is used to ex- 
plicate a connection between the infrared limit of the QCD Dirac operator's spectrum and the 
quark condensate appearing in the operator product expansion, and the connection is verified via 
comparison with a lattice-QCD simulation. The pseudoscalar vacuum polarisation provides a good 
approximation to the condensate over a larger range of current-quark masses. 
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I. INTRODUCTION 

Dynamical chiral symmetry breaking (DCSB) is a cor- 
nerstone of hadron physics. This phenomenon whereby, 
even in the absence of a current-quark mass, self- 
interactions generate a momentum-dependent running 
quark mass, M(p 2 ), that is large in the infrared: M(0) ~ 
0.5 GeV, but power-law suppressed in the ultraviolet 0: 
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is impossible in weakly interacting theories. In Eq. 
7 m = 12/(33 — 2Nf) is the mass anomalous dimension, 
with Nf the number of light-quark flavours, and (qq)° is 
the renormalisation-group-invariant vacuum quark con- 
densate , to which we shall hereafter refer as the OPE 
condensate. While Eq. JIJ is expressed in Landau gauge, 
(qq) is gauge parameter independent. In the chiral 
limit the OPE condensate plays a role analogous to that 
played by the renormalisation-group-invariant current- 
quark mass in the massive theory: it sets the scale of 
the mass function in the ultraviolet. 

The evolution of the dressed-quark mass-function in 
Eq. to a large and finite constituent-quark-like mass in 
the infrared, M(0) ~ 0.5 GeV, is a longstanding predic- 
tion of Dyson-Schwinger equation (DSE) studies Q that 
has recently been confirmed in simulations of quenched 
lattice-QCD 4]. A determination of the OPE condensate 
directly from lattice-QCD simulations must await an ac- 
curate chiral extrapolation [5j but DSE models tuned to 
reproduce modern lattice data give |(J (— {qq} ) ~ Aq CD . 

Another view of DCSB is obtained by considering the 
eigenvalues and eigenfunctions of the massless Euclidean 
Dirac operator 0: 



7 • Du n {x) = i\ n u n (x) 



(2) 



The operator is anti-Hermitian and hence the eigenfunc- 
tions form a complete set, and except for zero modes 



they occur in pairs: {u n (x) , j5U n (x)} , with eigenvalues 
of opposite sign. It follows that in an external gauge 
field, A, one can write the Green function for a massive 
propagating quark in the form 



S{x,y;A) = {q{x)q(y)) A 



E 



Un{x) u\(y) 
iX n + m 



(3) 



where m is the current-quark mass and, naturally, the 
eigenvalues depend on A. (NB. The expectation value de- 
notes a Grassmannian functional integral evaluated with 
a fixed gauge field configuration.) Assuming, e.g., a lat- 
tice regularisation, it follows that 
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where V is the lattice volume One may now define 
a quark condensate as the infinite volume limit of the 
average in Eq. |0J over all gauge field configurations: 



|m|0) := lim I f d'x ( {q{x)q(x)) A ) 

V Jy 



(5) 



In the infinite volume limit the operator spectrum be- 
comes dense and Eqs. (@J), ijjjj become 



-(0\qq\0) = 2m 



dX 



P(A) 



A 2 



(6) 



with p(X) the spectral density. This equation expresses 
an assumption that in QCD the full two-point massive- 
quark Schwinger function, when viewed as a function of 
the current-quark mass, has a spectral representation. 
It follows formally from Eq. © that 



hm ml dX -t. ^ 

m->0 J X z + m z 



~7TP(0), (7) 



and hence one arrives at the chiral limit result 

-(Q|gg|Q)=7rp(0). (8) 
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This is the so-called Banks-Casher relation ^J^J. It 
has long been advocated as a means by which a quark 
condensate may be measured in lattice-QCD simulations 
1 1 1| and has been used in analysing chiral symmetry 
restoration at nonzero temperature |12( and chemical po- 
tential |l3l ] , and to explore the connection between mag- 
netic monopoles and chiral symmetry breaking in U(l) 
gauge theory [lj| . Much has been learnt |15l Il6| by ex- 
ploiting the fact that qualitative features of the behaviour 
of p(A) for A ~ can be understood using chiral random 
matrix theory; i.e., from considerations based solely on 
QCD's global symmetries. 

Our main goal is to explicate a correspondence be- 
tween the condensate in Eq. Q and that in Eq. JSJ. 
In Sec. [H] we discuss the OPE condensate and its con- 
nection with QCD's gap equation, and emphasise that 
the residue of the lowest-mass pole-contribution to the 
flavour-nonsinglet pseudoscalar vacuum polarisation is a 
direct measure of the OPE condensate 17]. A natu- 
ral ability to express DCSB through the formation of 
a nonzero OPE condensate is fundamental to the suc- 
cess of DSE models of hadron phenomena ^^]. In Sec. 
IIIII we carefully define the trace of the massive dressed- 
quark propagator and use that to illustrate a connection 
between p(0) and the OPE condensate, which we verify 
via comparison with a lattice simulation. Section IIVI is 
an epilogue. 



II. OPE CONDENSATE 

A. Gap and Bethe-Salpeter equations 

Dynamical chiral symmetry breaking in QCD is readily 
explored using the DSE for the quark self-energy: 

S^)- 1 = Z 2 {i^ ■ p + m hm ) 

+ z 1 J^ ff 2 -D^(p-g)y7^(9)r^(g,p), (9) 

wherein: D^iji) is the renormalised dressed-gluon prop- 
agator; T^(q;p) is the renormalised dressed-quark-gluon 
vertex; m bm is the A-dependent current-quark bare mass 
that appears in the Lagrangian; and J := J d 4 g/ (27r) 4 
represents a translationally-invariant regularisation of 
the integral, with A the regularisation mass-scale which 
is removed to infinity as the completion of all calcula- 
tions. The quark-gluon-vertex and quark wave function 
renormalisation constants, Zi(( 2 , A 2 ) and Z 2 (C 2 , A 2 ) re- 
spectively, depend on the renormalisation point, the reg- 
ularisation mass-scale and the gauge parameter. 

If the current-quark mass changes with flavour, then 
the solution of Eq. © is flavour dependent: 

SAP)- 1 = i 1 -pA f {p 2 X 2 ) + B } {p\e) 

1 



Z f (p 2 ,C T 



[tj-p + Mfip 2 ^ 2 )] , (10) 



and is obtained subject to the condition that at some 
large, spacelike £ 2 

Sf(p)~\ 2 = C 2=il-p + m f (0, (11) 

where is the renormalised current-quark mass: 

Z 4 (C, A) m/(C) = Z a (C, A) m bm (A) , (12) 

with Z4 the renormalisation constant for the scalar part 
of the quark self-energy. Since QCD is an asymptotically 
free theory, the chiral limit is defined by 

Z 2 (C 2 ,A 2 )m^ m (A) =0, A>C (13) 

and in this case the scalar projection of Eq. © does not 
exhibit an ultraviolet divergence 0, El • 

Important in describing chiral symmetry is the axial- 
vector Ward-Takahashi identity: 
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-M (c) »Tf(fc;P)-ilf(fc;P)M (c)) 

(14) 

where P is the total momentum entering the vertex. 
In Eq. (HH): M (c) = diag[m„(C), m d (Q, m s (Q], S = 
diag[5 u , Sd, S 8 ] and {T H } are flavour matrices, e.g., 
T* =5 (A 1 + i\ 2 ) (we consider SU/(3) because chiral 
symmetry is unimportant for heavier quarks); r|^(fc;P) 
is the renormalised axial- vector vertex, which is obtained 
from the inhomogeneous Bethe-Salpeter equation 



[TUk;P)] tu = Z 2 



T 



H 



757m" 



(15) 



where X t ■= S(? + )rg,( ff ;P)S(g_), q± = q ± P/2, and 
K(q, k; P) is the fully renormalised quark-antiquark scat- 
tering kernel; and Tf? is the pseudoscalar vertex, 



[rf(* ; p)] t =z, 



T 



H 



75" 



[ X g(q;P)] ar K™(q,k;P), 



(16) 



with x" := S(q+)Tg(q;P)S(q-). Multiplicative rcnor- 
malisability ensures that no new renormalisation con- 
stants appear in Eqs. (fl5|l and l(Tojl. 

Flavour-octet pseudoscalar bound states appear as co- 
incident pole solutions of Eqs. I|15|) . I|16(l . namely, 



rf M (fc ; P) cx rf (fc; p) c< 1 g T H (k- P) 



(17) 



where Th is the bound state's Bethe-Salpeter amplitude 
and run, its mass. (Regular terms are overwhelmed at 
the pole.) Consequently, Eq. i|14|) entails |2ul7fl 



fiim H 



4°, 



(18) 
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where: M$ = tr flavour [M (c) {T ff , (T ff )*] is the sum of 
the constituents' current-quark masses ("t" denotes ma- 
trix transpose); and 



fnPu = Z 2 



r-A 



A 



1 



ir® = Z 4 J -tr 



1 



{T 11 )* TsT^feP)] , (19) 
(TBf-nxHtoP)] , (20) 



where xh S(q+)Tjj(q;P)S(q-) and the expressions 
are evaluated at P 2 + m 2 H = 0. 

Equation (|19fl is the pseudovector projection of the 
meson's Bethe-Salpeter wave function evaluated at the 
origin in configuration space. It is the precise expres- 
sion for the leptonic decay constant. The renormalisation 
constant, -^(C, A), ensures that the r.h.s. is independent 
of: the rcgularisation scale, A, which may therefore be 
removed to infinity; the renormalisation point; and the 
gauge parameter. Hence it is truly an observable. 

Equation (|20() is the pseudoscalar analogue. Therein 
the renormalisation constant Zi[(,, A) entails that the 
r.h.s. is independent of the regularisation scale, A, and 
the gauge parameter. It also ensures that the (- 
dependence of r^p is precisely that required to guarantee 
the r.h.s. of Eq. (|18|1 is independent of the renormalisa- 
tion point. (NB. r# is finite, and Eq. (|18(l valid, for 
arbitrary values of the current-quark masses jllj, l20| .) 

In the chiral limit the existence of a solution of Eq. 
@ with B (p 2 ) ^ 0; i.e., DCSB, necessarily entails [TtJ 
that Eqs. 1)15(1 . (|16[1 exhibit a massless pole solution: the 
Goldstone mode, which is described by 



Tg(fc;P) = A 9 7 5 



iE (k;P)+j-PF (k;P) 



(21) 



+ j-kk-PG (k;P) + a flv k^P l/ H (k;P) 



wherein f°E o (k;0) = B (k 2 ). (The index "0" indicates 
a quantity calculated in the chiral limit.) It follows im- 
mediately [TtJ that 



f o (0 

Jir '0 



■(qq)° c 



Urn Z 4 (C,A)iV c tr D / S (q) , (22) 

A ^°° Jo 



where the trace is only over Dirac indices. This result 
and multiplicative renormalisability entail 



<gg>g 



^4(C,C')^2" 1 (C,C')=^m(C,C') 5 



(23) 



where Z m is the mass-renormalisation constant. It is thus 
apparent that the chiral limit behaviour of yields the 
OPE condensate evolved to a renormalisation point C,. 

It is important to recall that the DSEs reproduce ev- 
ery diagram in perturbation theory. Therefore a weak 
coupling expansion of Eq. @ yields the perturbative se- 
ries for the dressed-quark propagator. This may be illus- 
trated by the result for the scalar piece of the propagator 



calculated in this way: 



B f (p 2 )=m f 1 



2 In 

7T 



(24) 



Every term in the series is proportional to the current- 
quark mass and hence a nonzero value of the OPE con- 
densate is impossible in perturbation theory. 



B. Pseudoscalar Vacuum Polarisation 

Consider the colour singlet Schwinger function describ- 
ing the pseudoscalar vacuum polarisation 

A 5 (z) = {q(x)±\ f l5 q(x) 9(0)^^759(0)} , (25) 

which can be estimated, e.g., in lattice simulations. Its 
renormalised form can completely be expressed in mo- 
mentum space using quantities introduced already: 

4 9 (P) = Z 2 tr J^T* %+) rf (g; P) S(g_) . (26) 

Equation (|16|) can be rewritten in terms of the fully- 
amputated quark- antiquark scattering amplitude: M = 
K + K(SS)K + . . ., and in the neighbourhood of the 
lowest mass pole 



M(q,k;P)=T H ( q] P) 



1 



P 2 



'H 



f(k;-P)+R(q,k;P), 
(27) 

where R is regular in this neighbourhood. 

Assuming SU/(3) flavour symmetry, substituting Eq. 
(ITl into Eq. J2gJ gives 



lj[ 9 (P) = 5 f9 



P 2 



7-2 2 



(28) 



(the ellipsis denotes terms regular in the pole's neigh- 
bourhood). It follows that the large- a; 2 behaviour of 

m 2 m A 5 (x) (29) 

is a measure of the renormalisation-group-invariant 



i(C) r£ 



(0 



(30) 



that appears in Eq. I|18jl. Hence the correlator in Eq. 
(|25|l provides a direct means of estimating the OPE con- 
densate in lattice simulations |2l|, one whose ultraviolet 
behaviour ensures a well-defined and calculable evolu- 
tion under the renormalisation group for any value of the 
current-quark mass. (NB. f„ can similarly be extracted 
from the axial- vector correlator analogous to Eq. JHJ).) 
The model of Ref. yields a meson mass trajectory via 
Eq. (|18fl that provides a qualitative and quantitative un- 
derstanding of recent quenched lattice simulations 20] . 
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III. BANKS-CASHER RELATION 



evaluated at a fixed value of the regularisation scale: 



A. Continuum analysis 

It is readily apparent that Eq. © is meaningless as 
written: dimensional counting reveals the r.h.s. has mass- 
dimension three and since A will at some point be greater 
than any relevant internal scale, the integral must diverge 
as A 2 , where A is the regularising mass-scale. 

To learn more, consider the trace of the unrenormalised 
massive dressed-quark propagator: 



cr(m) := N c tro / S m (p) , 
Jp 



(31) 



evaluated at a fixed value of the regularisation scale, A. 
This Schwinger function can be identified with the l.h.s. 
of Eq. |JBJ. Furthermore, assume that a(m) has a spectral 
representation, since this is the essence of the Banks- 
Casher relation: 



a(m) := 2m f d\- ^ 



A 2 



where m = m hm (A). Equation ilML'l) entails 



(32) 



p(A) = — lim [ a(iX + rj) - a{i\ - V ) ] . (33) 

The content and meaning of this sequence of equations 
is well illustrated by inserting the free quark propagator 
in Eq. I|31|l . The integral thus obtained is readily evalu- 
ated using dimensional regularisation: 



r n N c 3 



In — — H h 7 — m 4ir 

C £ 



(34) 



With Eq. 1331) the regularisation dependent terms cancel 
and one obtains 



P(A) 



4tt 2 



A 



(35) 



The one-loop contribution to /5(A) has been evaluated 
using the same procedure |23j |. It is also proportional 
to A 3 and arises from the m 3 lnm 2 /( 2 terms in a(m). 
In fact, every term obtainable in perturbation theory is 
proportional to A 3 , for precisely the same reason that 
each term in the perturbative expression for the scalar 
part of the quark propagator is proportional to m, see 
Eq. i|24|) . Hence, at every order in perturbation theory, 



p(A = 0) = 



(36) 



and (0\qq\0) — 0. A nonzero value of p(0) is plainly an 
essentially nonperturbative effect. 

A precise analysis requires that attention be paid to 
renormalisation. Consider then the gauge-parameter- 
independent trace of the renormalised quark propagator 



cr(m;C) := ~ \im (q(x)q(0)) 

x—*0 



= Z 4 (C,A)iV c tr D / S m (p;0, (37) 
Jp 

where the argument remains m — m bm (A), which 
is permitted because m bm (A) is proportional to the 
renormalisation-point-independent current-quark mass. 
The renormalisation constant Z± vanishes logarithmi- 
cally with increasing A and hence one still has a(m) ~ 
A 2 m bm (A). However, using Eq. (|2*2"|) it is clear that for 
any finite but large value of A and tolerance 5, it is always 
possible to find m^ m (A) such that 



a(m; C) + (qq)° c < S , Vm bm < m bm . 



(38) 

This is true in QCD. It can be illustrated using the 
DSE model of Ref. 0, which preserves the one-loop 
renormalisation group properties of QCD. In Fig. ^ we 
plot cr(m, C), evaluated using a hard cutoff, A, on the 
integral in Eq. H37[). calculated with the massive dressed- 
quark propagators obtained by solving the gap equation 
as described in the appendix. Since Eq. specifies the 
domain on which the value of <r(m; £) is determined by 
nonperturbative effects, one anticipates 



bm 



(A^-^/A^lO- 



for A = 2.0 TeV in QCD where \{qq)°\ 
mate confirmed in Fig. ^ 
The dotted line in Fig. ^ is 

2 A 
cr(m, Q = — (qq)® — arctan — 



A 3 



(39) 
an esti- 



[A 2 -m 2 ln(l + A 2 /m 2 )] .(40) 



(NB. We used the one-loop formula: Z^CA) = 
[a(A)/a(C)] 7 "*, for the numerical comparison.) The 
difference between Eq. (|40ll and the curve is of 
0(a(A) m(A) A 2 ) because the DSE model incorporates 
QCD's one-loop behaviour. In Fig.^we also plot a(m, () 
obtained in the absence of confinement, in which case [24| 
(qq)° = 0, as is apparent. 

The discussion establishes that a(m, £) has a regular 
chiral limit in QCD and is a monotonically increasing 
convex-up function. It follows that <r(m, () has a spectral 
representation: 



= 2m f 
Jo 



A 2 + m 2 



(41) 



This lays the vital plank in a veracious connection be- 
tween the condensates in Eqs. and (jHJ. On the do- 
main specified by Eq. I|39|l . the behaviour of er(m, £) in 
Eq. (J23 is given by Eq. gUJ), which yields, via Eq. (|33|> . 



Trp(A) 



-(m)° c 



Z 4 (C,A)^A 3 + 



(42) 



where the ellipsis denotes contributions from the higher- 
order terms implicit in Eq. I|40|) . 
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FIG. 1: Circles/solid-line: cr(m) 1,/3 in Eq. I|M7|I as a function 
of the current-quark bare-mass, evaluated using the dressed- 
quark propagator obtained in the model of Ref. |3; dashed 
line: the model's value of (— (qq}° =1 Ge v) = (0.24 GeV) 3 ; and 
dotted line: Eg. 14UII . Diamonds: a(m) 1,/3 evaluated in a non- 
confining version of the model; dot-dashed line: Eq. 14011 with 
(qq)° = 0. 
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FIG. 2: Spectral density of the staggered Dirac operator 
in quenched SU(3) gauge theory calculated on a 4 4 -latticc. 
(Measured in units of the lattice spacing.) The deconfinement 
transition takes place at (5 > 5.6. 



B. Comparison with a lattice-QCD simulation 

In Fig. [21 we plot the spectral density of the staggered 
Dirac operator in quenched SU(3) gauge theory calcu- 
lated with 3000 configurations obtained on a V = 4 4 - 
lattice, in the vicinity of the deconfining phase transition 
at 8 > 5.6. Details of the simulation are given in Ref. 
[12) . Dimensioned quantities are measured in units of 
1/a, where a is the lattice spacing, and it is p(X)/V that 
should be compared with the continuum spectral density. 

While the effect of finite lattice volume is apparent in 
Fig- for Xa > 0.1, the behaviour at small Xa is qualita- 
tively in agreement with Eq. I|42|l and Fig. a nonzero 
OPE condensate dominates the Dirac spectrum in the 
confined domain; and it vanishes in the deconfined do- 
main whereupon p(0) = and the perturbative evolu- 
tion, Eq. I|35[) . is manifest. 

To be more quantitative, we note that at = 5.4, 
p(0)a 3 w 70, so that 

np{0)a 3 /V w (0.95) 3 . (43) 

The value of the lattice spacing was not measured in Ref. 

but one can nevertheless assess the scale of Eq. (|43f) 
by supposing a ~ 0.3 fm ~ 0.3/Aqcd, a value typical of 
small couplings, f3, wherewith the r.h.s. is ~ (3Aqcd) 3 - 
This is too large but not unreasonable given the parame- 
ters of the simulation, its errors and the systematic uncer- 
tainties in our estimate. One can also fit the lattice data 
at (3 — 5.8, whereby one finds p(X) oc A 3 on A < 0.1 but 
with a proportionality constant larger than that antici- 
pated from perturbation theory; viz. Eq. I|42|) . Some mis- 
match is to be expected because at (3 = 5.8 one has only 



just entered the deconfined domain and close to the tran- 
sition boundary nonperturbative effects are still material, 
as seen, e.g ., in the heavy-quark potential and equation 
of state [2{| • It is a modern challenge to determine those 
gauge couplings and lattice parameters for which the data 
arc quantitatively consistent with Eq. I|42|l . 

IV. EPILOGUE 

We verified that the gauge-invariant trace of the mas- 
sive dressed-quark propagator possesses a spectral repre- 
sentation when considered as a function of the current- 
quark mass. This is key to establishing that the OPE 
condensate, which sets the ultraviolet scale for the 
momentum-dependence of the trace of the dressed-quark 
propagator, does indeed measure the density of far- 
infrared eigenvalues of the gauge-averaged massless Dirac 
operator, a la the Banks-Casher relation. This relation 
is intuitively appealing because a measurable accumu- 
lation of eigenvalues of the massless Dirac operator at 
zero-virtuality expresses a mass gap in its spectrum. 

In practice, there are three main parameters in a sim- 
ulation of lattice-QCD: the lattice volume, characterised 
by a length L; the lattice spacing, a; and the current- 
quark mass, m. So long as the lattice size is large 
compared with the current-quark's Compton wavelength; 
viz., L 3> 1/tj, then dynamical chiral symmetry breaking 
can be expressed in the simulation. Supposing that to be 
the case then, as we have explicated, so long as the lat- 
tice spacing is small compared with the current-quark's 
Compton wavelength; i.e., a <§; l/m <ti L, 

np(0)K-(qq)1 /a , (44) 
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where the r.h.s. is the scale-dependent OPE condensate 
[C = A = l/o in Eq. g2J]. 

In our continuum analysis we found that one requires 
am < ((zAqcd) 3 if p(A = 0) is to provide a veracious esti- 
mate of the OPE condensate. The residue at the lowest- 
mass pole in the flavour-nonsinglet pseudoscalar vacuum 
polarisation provides a measure of the OPE condensate 
that is accurate for larger current-quark masses. 
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APPENDIX A: MODEL GAP EQUATION 

The gap equation's kernel is built from a product of the 
dressed-gluon propagator and dressed-quark-gluon ver- 
tex. It can be calculated in perturbation theory but that 
is inadequate for the study of intrinsically nonperturba- 
tive phenomena. To make model-independent statements 
about DCSB one must employ an alternative systematic 
and chiral symmetry preserving truncation scheme. 

The leading order term in one such scheme pr| is 
the renormalisation-group-improved rainbow truncation 
of the gap equation (Q = p — q): 



S^y 1 = Z 2 (i 1 -p + m hln ) 



G{Q z W™ c {Q)^i»S{ q )- lv 



(Ai) 



The ultraviolet (Q 2 > lGeV 2 ) behaviour of £(Q 2 ) in 
Eq. (|A1() is fixed by the known behaviour of the quark- 
antiquark scattering kernel Q ■ The form of that kernel 
on the infrared domain is currently unknown and a model 
is employed to complete the specification of the kernel. 
An efficacious form is Q 



= 8ir 4 D5 4 (k) 



Q 2 

+47T- 



-DQV 



7m 7T 



■hi 



? 2 /A 2 



QCD 



HQ 2 ), (A2) 



where: T{Q 2 ) 
0.5 GeV; r = e 2 

A A W _ 



1 



[1 

7r> 



exp(-Q 2 /[4m 2 ])]/Q 2 , m t = 
12/(33 -2Nf), N f = 4; and 
qcd = A^i = 0.234 GeV. The true parameters in Eq. 
(|A2|) are D and oj, however, they are not independent: 
in fitting, a change in one is compensated by altering 
the other, with fitted observables changing little along a 
trajectory ujD — (0.6 GeV) 3 . Herein we used 



D = (0.884 GeV) 2 



0.3 GeV. 



(A3) 



A non-confining model is obtained with D = 0. 

Equation (|Alj) . is readily solved for the dressed-quark 
propagator, with the renormalisation constants fixed via 
Eq. I|ll|) ■ That solution provides the elements used in the 
illustration of Sec. IIIII 
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